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ABSTRACT 

A new  approach  is  supplied  for  evaluating  tin*  efficiency 
of  decision  making  units,  locating  efficiency  frontiers  and  estimating 
parameter*  from  observational  data.  This  la  accomplished  by  meana  of 
a nonllnear-nonconvex  hl-extremal  principle  which  In  subsequently 
shown  to  be  essentially  reducible  to  a finite  sequence  of  linear  programming 
problems.  The  development  la  Illustrated  by  means  of  multiple  output 
functions  which  are  piecewise  of  Cohb-Douglas  or  general  log  linear  type 
and  which  also  allow  for  Increasing,  decreasing  and  constant  returns  to 
scale.  The  reduction  of  the  bi-extremal  principle  to  linear  programming 
equivalence  la  also  accomplished  for  much  more  general  classes  of  functions. 


A 
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1.  INTRODUCTION 

CCR  13)  provided  a nonlinear  ratio  extremal  principle  for  determining 
a piecewise  linear  efficiency  frontier  for  a collection  of  DMU's  (Decision 
Making  Units).  This  is  fundamental  to  the  subsequent ly  developed  procedures 
of  Data  Envelopment  Analysis  which  are  used  to  distinguish  between  "Program," 
"Managerial"  and  other  types  of  efficiency.  See  [2).  One  limitation 
of  the  preceding  analysis,  however,  is  that  the  efficiency 

frontiers  are  formed  from  functions  which  are  of  linear  and/or  piecewise 
linear  type.  Some  situations  may  Involve  other, more  general, parametric 
families  of  functions.  One  may  then  need  to  estimate  the  parameters 
with  assurance  that  the  resulting  functions  give  the  relevant  efficiency 
t ront iers. 

Here  we  shall  develop  a bi-extremal  principle  for  simultaneously 
achieving  the  parametric  estimates  and  the  associated  efficiency  measure  to 
be  assigned  to  each  DMU.  The  basic  idea  is  as  follows:  We  "envelope" 
the  observed  values  of  the  outputs  by  means  of  parametric  functions  of 
the  observed  input  values.  The  envelope  is  "tightened"  to  rest  on  possible 
observed  output  values  by  means  of  a minimization  operation.  Then  the 
parametric  output  functions  are  used  to  replace  the  observed  output  values 
in  the  manner  of  the  nonlinear  ratio  maximization  principle  of  CCR  (3) 
tor  the  efficiency  determination  of  each  DMU.  Thus  one  obtains  a bi-extremal 
principle  of  "mpximim  type." 

In  this  paper,  and  for  output  functions  of  Cobb-Douglas  type 
(or,  more  generally,  for  "<J>  linear"  type)^  we  show  that  this  nonlinear- 
nnneonvex  bi-extremal  principle  is,  in  fact,  equivalent  to  a finite  sequence 
of  linear  programming  problems  with  one  nonlinear  vector  operation.  The 
linear  programs  further  differ  only  in  the  right-hand  side,  thus  greatly 
facilitating  computation  and  analysis. 

1]  See  sect  ion  6 below. 
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Ko»t  current  approaches  to  the  estimation  of  extreme  values 
(a.g.,  as  In  mathematical  statistics)  are  restricted  to  *ero  dimensional 
values  (or  frontiers)  such  as  the  end  point!  of  a uniform  distribution. 

Our  bl-extremal  principle  provides  a new  approach  to  the  determination 
of  multi-dimensional  (extreme)  frontiers  such  as  are  required  for  efficiency 
determination  in  the  multiple  Input-multiple  output  situations  that 
are  comswnly  encountered  in  public  policy  evaluation  problems.  The 
principle  la  more  general  than  such  aniltlple  input-multiple  output 
applications,  however,  and  extends  to  any  case  where  extreme  frontiers 
are  of  Interest. 

2.  BACKGROUND 

Wc  follow  the  notation  conventions  of  13)  and  consider  the  coamom 
Input  and  output  values  for  a collection  of  DMU’s  defined  as  follows: 


(1) 


x^  • the  amount  of  Input*  l • l,  m; 

y^  " the  aoxHint  of  output,  r ■ 1,  ....  s. 


where  ) • 1 n Indexes  each  one  of  the  DMU*s  betn*  considered. 

These  x^  and  values  will  generally  represent  observations  generated 

from  past  behavior  and  we  shall  assume  that  they  all  have  positive  values. 

The  following  formulation  was  given  In  |3]  for  determining 
the  efficiency  of  any  specified  DMUq  from  among  this  set  of  J - l DMU's 
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Because  the  ratio  In  the  functional  also  appears  In  the  constraints  we 

have  mas.  h • h*  ^ 1 In  any  case  ami,  as  shown  In  |3|,  h*  ■ 1,  Iff 
0 O’"  o 

UMV)^  Is  efficient.  Note,  In  particular,  that  this  provides  a scalar  value 
for  the  wanted  efflciencv  rating. 

Thla  scalar lxat Ion  Is  achieved  via  the  non-negative  weights  w^, 
assigned  to  the  respective  outputs  and  inputs.  These  weights  are  not 
assigned  g mloil  In  an  arbitrary  manner,  however,  but  are  determined 
objectively  from  the  data  as  prescribed  by  (2).  Hereafter 

we  shall  refer  to  them  as  "virtual  weights  " — the  intended  analogy  being 
to  concepts  like  "virtual  displacements"  and/or  "virtual  work"  In  physics 

or  engineering  , which  represent  magnitudes  that  are  not  observed  directly 
but  ate  Implicit,  Instead,  tn  the  underlying  pnyslcal  principles  a no  models. 

In  the  present  paper  we  want  to  replace  the  formulation  In  (2) 
with  another  more  general  one  that  will. inter  alia,  enable  us  to 
estimate  efficiency  for  t'obb-Pouglas  type  (multiple  output)  functions.  We 
want  to  do  this,  however,  without  losing  contact  with  the  developments  in  13). 
In  part tcular, we  want  tv'  preserve  the  Data  Envelopment  Analysis  procedures 
of  the  preceding  work  so  that  we  will  thereby  be  able  to  effect  evaluations 
that  distinguish  between  programs  or  technologies  and  managerial  efficiencies 
that  are  of  importance  for  public  sector  guidance  and  control.  See  (2). 

To  help  fix  the  Ideas,  we  present  our  development  In  terms  of  Cobb-Douglas 
functions  and  then  present  the  more  general  formulae. 


Jl, 


'A  discussion  of  the  operational  significance  of  this  rating  as  well 
as  a t ransformat Ion  into  linear  programming  equivalents  (o.g.  lor  computational 
efficiency)  are  provided  In  [ ). 

Jj 

See,  e.g.,  the  discussion  on  p.b47  in  (1). 


3.  DEVELOPMENT 


First,  we  Introduce  th*  " envelopment  condition"  on  the  outputs. 


”rj  i \ ,2,  «,J  ' 


.1  * ^ • i n i 


where  the  variable-  and  Af  sre  const  rained  to  be  non-negative. ^ Here 

* u 

the  symbol  ^ reteis  to  the  *-fold  product  of  the  rt  vsluea  so 

that  we  are  now  restricting  attention  to  functions  which  ore  of  Cobl>- 


Douglas  form  — but  extended  to  the  case  of  multiple  outputs  t • 1 s 

In  number . 


A*  already  observed , we  wish  to  make  these  iih  qualities  aa  tight 
as  possible.  Hence  we  orient  our  ob|ectlve  toward  a minimisation  over 


the  A^  and  variables,  as  In  the  following 


max.  min. 


s m U 

E.w  A n.x/1 

r-1  r r r-1  lo 
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m it 

y < A n.x.T1 
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Recall  that  y and  x are  observed  positive  constants. 


See  1 4 J and  (5)  lot  further  extension  to  the  class  of  positive  homogeneous 
and  analytic  functions  for  these  Cohh-Douglas  tvpe  of  formulat Ions. 


litis  formulation  provides  our  non  1 Ineir-noiicoitvex  lil-txtrem.il  |’iincl|>le. 
here  formulated  for  functions  oi  Cobb-Douglna  (multiple  output)  type.  Note 

that,  unlike  the  development  in  |»1,  we  here  have  a a i mu 1 1 ancons  determination 
°*  (1)  the  parameter  values  A^,  ^ and  (11)  the  efficiency  rating"  via  (ill) 

the  "virtual  weights"  w and  u . 

r * 

The  formulation  in  (4)  involves  non- . ouv.  s U v (or  non-c.  nenvityi 

...  the  bi-extremal  objective.  However,  a-,  not . d ...  Introduction. 

we  shall  bring  it  into  an  equivalent  concave-convex  torn,  reducible  to 
an  a_l’r±orl  I ixed  I luite  sequence  ol  linear  piogi  miming  problems  with 
one  nonlinear  veetei  operation. 


4 . TRANS  IMKMAT I ON 


We  new  introduce  the  variables  z via 

ro 


O) 


r i 


ro 


r 1»1  .o  - 

As  a result  of  adjoining  these  additional  const. dints  to  (4)  an  equivalent 
objective  function  is 


(«’) 
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we  make 

the  following  change  oi  vai iahle 
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• ( w 

r 
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such  that 
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In  those  now  variables  wo  have  .1  blltno.ir  functional  and  vOnatralnts 


••  follows 


MX.  win.  r w'.* 
• • r*l  r 
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By  taking  logs  in  the  first  two  lna^oal it io#  wo  onn  rewrite  (9) 


in  the  form. 


max.  min.  l.w'z 


r-l  r ro 


ro 


subject  to 


(10) 


m 

r 

i-1 

“rl 

*lo 

^ 0 

m 

i-ri 

li  . 

r 1 

X11 

< -V 

i*lul  *tJ 


— I.  w v . - .I.U.  X., 

r*l  r rj  i“l  l ij 


1 

0 


"ro*  Ar*  Mrl.  v^.  u^  > 0 


where  the  caret  above  a symbol  denotes  natural  logarithm. 
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*>.  LINEAR  PR0C.RAMM1NC  EtflJ  1VAI.F.NT 

Notice  In  (10)  that  the  w',  u'  and  the  * , A , u sets  of 

ro  r rl 

variables  occur  In  scparate,lndependent  seta  of  constraint*.  These  two 

constraint  systems  restrict  their  separate  v.u  table  seta  to  convex  and 

feasible  sets.  Further,  with  the  bilinear  functional,  the  bl-extremal 

problea  la  equivalent  to  a saddle-value  problem  and  such  an  equivalence 

would  hold  for  much  more  general  bl-extremal  problems  than  (10). 

We  shall  pursue  this  generality  in  subsequent  papers. 

Here  we  note  next  that  the  ? , A , p , constraint  system 

ro  r rl 

separates  Into  independent  systems  with  one  system  tor  each  output  Index  r. 

Thus  for  any  set  ot  the  non-neg.it  ivo  w^  values,  we  wish  to  minimize  each 

* _ subject  to  its  separate  (.  r * *'  > constraint  system.  Setting  z - fn(z  ), 
ro  " ro  ro 

l he  r problem  would  be: 

* A 

X 

. I O 

min.  e 


ild.l)  subject  to 


r ~ A - .u.H  .x,  > 0 

ro  r ri  lo  - 


A ♦ — it  \ v , |"1 n 


i*r  ri  I | 


r ) 


“ri  0 ’ 


i I . ....  m. 


Or,  since  mlniratz  tup.  e is  equivalent  to  mintmi.’ing  ? , we  have  the  lin  ar 

ro 

pt  k g ramming  problem. 


m i n . 


ro 


ill) 


subject  to 
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« i lHi  i I 1 


i V 


i 0 • m . 


Notice  that  the  functional  and  matrix  coelf lclent  structure  of 


this  linear  programming  system  Is  the  same  for  all  r-1.  ....  s.  Only  the 
right  hand  side,  the  y ^ , would  change  with  r.  II  one  were  to  solve  the 
system  via  the  dual  linear  programming  problems,  only  the  functional 
would  change  with  r.  One  then  knows  a priori  each  functional  to  be 
employed.  Thus  an  optimal  has  >;  for  one  value  of  r would  he  a feasible 

, we  make  the  inverse 
thereby  secure  the  following 
, u'*,  (hence  , u*  if  desired): 

- 1 


nas is  i or  any  vam*  01  r 


On  obtaining  the  optimal  z rt* 

ro  ro 

(exponential)  t ransform.it  Ion  to  get  z*  and 

ro 

linear  programming  problem  to  determine  w'* 


i.ix.  . w’  z* 
r 1 r ro 


(12) 
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, • , u ' x 
1*1  i lo 


subject  to 


r*  1 UrVr J 


— J U * X 

l- l 1 i l 
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The  mjxlmal  value  of  the  functional!  :.w*'z*  )is  then  t lie  efllciencv  rating 

r ro 

for  llMl’  . .similarly,  the  optimal  values  of  the  variables, 

m lirl 

Ujj  give  us  tlu  local  Oobb-Douglas  type  envelope  A*  function 

at  this  »otn"  DMU. 


4 cenkhai.i/ation 


It  is  clear  that  to  get  equivalent  linear  programming  systems  of  the 

sort  of  (11)  and  (12)  to  our  bl-extren*al  principles  much  more  general  output 

functions  than  those  of  Cobb-lhuiglas  or  genet. il  log  linear  tvpe  are  permissible. 

We  shall  call  an  output  f one  t Ion  • . . ,#  , x ....  x ) "♦-linear"  Iff: 

1 *■ 

ml 

<n>  <i|g(t,x)|  * eKfl^x) 

where  ( Is  a mono! me  strictly  increasing  function. 

We  shall  assume  each  output  function  g^  is  j^-llnear.  Applying 

the  same  techniques  as  In  our  Cohh-|)ouglas  cases,  we  can  reduce  our  bl-extrss»sl 

principle  to  solution  of  the  following  w,... 

toiiowfnjt  two  linear  programming  systems.  First, 


min. 


n. 

'•  , f ,l»  ) 

ro  k*l  rk  rk 


04. 1) 


Dufijvct  to  rk*  rk^x  ^ — *r^v  ^ ^ • ••••  n» 


pin*  ,uii  other  relevant  linear  iiuiiii.i  1 1 1 ies  on  (In  9 wliiili  do  not 

rk 

involve  am  other  value  ol  r.  I'Mi  min  . •'.*  available,  the  efficiency  system 

ro  ro 

t hen  becomes  tin  s.  mnl  ol  these  two  linen  |>rmi  mis  - v i / . , *■ 


max.  . w ' ; * ( * ) 
r«l  r ro 


(14.2) 


subject  to 


in 

i i 

v*  y . • £u.  n 
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f * 

w , u > 0 • 
r i — 


CONCLUSION 

Between  CCK  [3]  and  the  above, we  have  now  explicitly  covered  the 

ease  ol  ^-linear  functions  — winch  includes  log  linear  and  Cobb-Douglas 

families  ot  functions.  Convenience  of  use  and  the  extent  of  their  validation 

in  empirical  studies  makes  members  of  our  Cobb-Doug : as  family  of  special  interest.lt 

should  be  noted,  moreover,  that  our  formulations  extend  to  multiple  outputs  and 

to  iiiic'ency  frontiers  wnicn  are  piecewise  Cobb-Douglas,  Also, 

we  did  not  rest  riit  the  choice  ol  u .to  the  cast,  of  constant  returns  to 

' r i 

scale.  Hence,  unlike  the  piecewise  linear  case  of  | 1 1 , we  now  permit 
in  n.ising  and/or  decieasing  returns  to  scale  both  in  the  frontiers  and 
in  the  PMC's  being  evaluated. 

The  measures  for  efficiency,  howevi r,  retain  their  operational 
significance.  See  [ J|.  It  should  he  spec i f ica ' l v noted  that  in  the  eases  of 
increasing,  decreasing  and  constant  returns  to  scale  the  etficienev  evaluation 
is  effected  by  reference  to  the  frontier  segments  that  have  this  same  property.  In 


other  words,  we  can  now  ascertain  whether  a DMU  that  has  Increasing 
returns  to  scale  is  operating  efficiently  within  any  such  zone  as  well 
as  whether  it  is  attaining  scale  economies.  Moreover,  we  can  also 
restrict  particular  coefficients  so  that  they  provide  only  one  of 
Increasing,  decreasing  or  constant  returns  to  scale.  This  could  be  done 
by  lntroducting  inequalities  on  the  exponents  or  functions  without 
altering  any  of  the  properties  of  the  above  models.  Finally  we  should 
note  that  these  multiple  output/multiple  input  formulations  permit 
these  features  to  be  used  in  some  of  the  output/input  relations  without 
requiring  them  for  all  of  the  others  at  the  same  time. 
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